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In ferromagnets, charge currents can excite magnons via the spin-orbit couphng. We de- 
velop a novel and general scattering theory of charge current induced macrospin magnetization 
torques in normal metal |ferromagnet| normal metal layers. We apply the formalism to a dirty 
GaAs|(Ga,Mn)As|GaAs system. By computing the charge current induced magnetization torques 
and solving the Landau-Lifshitz-Gilbert equation, we find magnetization switching for current den- 
sities as low as 5 X 10^ A/cm^. Our results are in agreement with a recent experimental observation 
of charge- current induced magnetization switching in (Ga,Mn)As. 



Magnetic random access memories are read using the 
principle of the giant magnetoresistance effect. Informa- 
tion is written by non-local magnetic fields, and this lim- 
its the density of the memory cells. In small systems, a 
promising direct and local writing technique is to utilize 
spin-transfer torques from spin-polarized currents in non- 
collinear magnetic structures [1 . The interplay between 
out-of-equilibrium quasi-particle spin-polarized currents 
and collective magnetization dynamics exhibits a rich va- 
riety of phenomena, and receives currently much atten- 
tion from the scientific community. When a spin po- 
larized current enters a ferromagnet, the itinerant flow 
of angular momentum normal to the magnetization di- 
rection is absorbed. The absorbed angular momentum 
is transferred to the ferromagnetic order parameter as a 
torque on the magnetization [1 . Magnetization switch- 
ing has been seen in several experiments at current den- 
sities of the order 10^ - 10^ A/cm^ [1 . 

Most theories of spin-transfer torques in layered sys- 
tems disregard spin-orbit interaction in the band struc- 
ture and impurity potentials. Interestingly, the spin- 
orbit interaction, which is responsible for spin memory 
loss, also enables transfer of orbital angular momentum 
to the collective magnetization. Thus, an unpolarized 
current can induce a torque on the magnetization even 
in magnetic homogeneous systems. This has recently 
been demonstrated theoretically for a simplified mag- 
netic 2DEG Rashba Hamiltonian [2 and experimentally 
seen in (Ga,Mn)As [3 . The experimental results call 
for 1) an extension of current induced torque theories to 
include spin-orbit interaction which is valid for general 
band-structures and disorder as well as a 2) quantita- 
tive computation of the current induced magnetization 
dynamics. 

We consider a ferromagnet in electric contact with nor- 
mal metals. The torque on a macrospin ferromagnet 
characterized by a spatially independent unit vector m 
can be written as 

m = -7m X [Heff(m) + H,(m) + Hc(m)] , (1) 

where 7 is the gyromagnetic ratio and Heff is the effective 



magnetic field arising from the ferromagnet 's equilibrium 
energy. H5 and He are induced by non-equilibrium quasi- 
particle transport. Hg oc 11 ^ describes the spin-transfer 
torque induced by a spin accumulation in the normal 
metals adjacent to the ferromagnet. Our main focus is 
the charge-current induced torque, represented by He cx 
induced by a voltage bias V across the ferromagnet 
without any spin bias. When the spin-orbit interaction 
vanishes so does He. 

We will in this paper derive expressions for the field He 
valid for general bandstructures and disorder, and with 
that formulate a novel and general theory for current- 
induced magnetization dynamics in layered N|F|N struc- 
tures. The theory describes how unpolarized charge cur- 
rents can excite magnons in the presence of spin-orbit 
coupling. We also find generalized expressions for Hg 
that is valid even when there is spin-orbit coupling and 
the spin-transfer torques cannot just be evaluated from 
the loss of transverse spin current. Without spin-orbit 
coupling. He = and our formalism reduces to the mag- 
netoelectronic circuit theory [4 . We apply our theory 
to a disordered GaAs|(Ga,Mn)As|GaAs system, and find 
that an unpolarized charge current induces magnetiza- 
tion reversal for current densities as low as 5 x 10^ A/cm^. 
This is in agreement with recent experimental results and 
demonstrates the usefulness and potential of the scatter- 
ing theory. 

Ferromagnetic semiconductors are examples of systems 
with a large intrinsic spin-orbit coupling which opens new 
paths to the manipulation of the magnetization direction. 
In (Ga,Mn)As, ferromagnetism is mediated by holes, and 
the magnetic anisotropy energy depends on the hole con- 
centration. This together with the ability to control the 
hole density by gate voltages in semiconductors, enables 
magnetization switching [5 . In contrast, in our case, it 
is the current induced out-of-equilibrium hole population 
that initiates magnetization switching. 

Spin is not conserved in systems with spin-orbit cou- 
pling so spin transfer torques cannot be found from the 
discontinuity of the transverse spin-current at normal 
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FIG. 1: (Color online) (a) and (b) are reciprocal spin pro- 
cesses: (a) spin transfer torque induced magnetization dy- 
namics by a spin bias /itf (m? ) ^he left (right) reservoir, 

(b) spin pumping by a precessing magnetization into the left 
and right reservoir, (c)-(d) are reciprocal charge processes, 

(c) magnetization dynamics induced by voltages Vl and Vr, 

(d) charge pumping by precessing magnetization between the 
left and right reservoir. Observation of one of the processes 
implies, by Onsager reciprocal relations, the existence of the 
reciprocal process. 



metal|ferromagnet interfaces. Instead, torques can be 
formulated in terms of the spatially dependent exchange- 
correlation energy [6 , but this requires a computation of 
the out-of-equilibrium spin density in the ferromagnet. 
We choose another route using the Onsager reciprocal 
relations resulting in a compact result for the torques in 
terms of the scattering matrix of the complete system 
without requiring a detailed knowledge of the possibly 
complicated spatially dependent out-of-equilibrium spin 
density. 

The derivation of our formalism is based on the 
Onsager reciprocity relations, which state the follow- 
ing \T/. Assume the system is described by the param- 
eters {qi\i = 1, . . . , N}, and the rate of change Qi is in- 
duced by the thermodynamic force Xi . In linear response 

Qi = ^j^=i^ij^j- ^^^^ ^ij = ^i^jLji, where 6^ = 1 
(ci = —1) if Qi is even (odd) under time reversal In a 
N|F|N system a voltage V^(^) and a spin accumulation 
/jLs^^"^ at the left (right) node induce a rate of change 
in the charge and spin density in the left (right) reser- 
voir, which by the continuity relation equals the charge 
and spin flows and if^^^ into the ferromagnet. 

Xc = {V^ — V^) y is the force on the charge system 
inducing the charge current I = /y along the transport 
direction y as illustrated in Fig. [TJc). Similarly, the spin 
force xf^^'^ = ^^^^^ /% at the left (right) node induces 
the spin current Fig. Il[a). The magnetic system 



is described by the magnetic free energy F[M] and the 
macrospin M = l/Mgin, where V is the volume of the 
ferromagnet, Ms the absolute value of the magnetization, 
and m the unit direction vector. In this case, the effec- 
tive field Xm = — 9m^[M] induces the response M in 
the magnetization. The relations between the rates and 
the thermodynamic forces are 
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where L^^^ ^ characterizes the spin transfer torque arising 
from a spin accumulation in the left (right) node, Lmc 
represents the torque from a charge potential, Lmm de- 
scribes the magnetization dynamics for an isolated mag- 
netic system according to the Landau-Lifshitz equation. 

We identify the relations [Lms = [^sM 

and [Lmc (m)] - [Lcm {-m)]... The (Lsm)^^- matrix 
describes spin pumping by a precessing magnetization. 
Its reciprocal process is spin transfer torque induced by 
a spin accumulation, as illustrated in Fig. [TJa)-(b). Spin 
pumping is a well known phenomenon predicted theoreti- 
cally [8 and measured in experiments [9 . Charge pump- 
ing is described by (Lcm)^^ and caused by the spin-orbit 
interaction. Charge pumping is the reciprocal process of 
charge potential induced magnetization dynamics, and is 
therefore a reciprocal manifestation of a charge current 
torque and non- vanishing coefficients (LMcJij- This is 
illustrated in Fig. [T]^c)-(d). 

The coefficients L^^'' and Lmc can be found from the 
out-of-equilibrium spin density in the ferromagnet, but 
such expressions are lengthy and cumbersome when spin- 
orbit coupling is taken into account. Instead, we consider 
the reciprocal processes, charge and spin pumping by 
magnetization precession, and compute the coefficients 
L^l^^ and Lcm- Parametric pumping can be found from 
the scattering matrix of the N|F|N system ^ and using 
Mi {LMM)ik (Xm)/c it reads: 

i^(^)=Qf^^\LMMhA^M)k. (3) 

Qk = y Yl ^mn,wdmjlnn,W + ^.C. , (4) 

where Smn,w is the 2x2 scattering matrix in spin- 
1/2 space from transverse mode n in lead V to trans- 
verse mode m in lead The current operator /^(^) 
is in the 2x2 spin space which can be expanded in 
terms of charge and spin current contributions, = 
^q/^(^)/2 + ea • Is^^"^ /h^ where (Jq is the unit ma- 
trix and & is di vector of Pauli matrices. This rep- 
resentation and Eq. ([s]) give us the response coeffi- 
cients {LcM)yk = Tr[Qf - Q^]{Lmm),J2 [TOl and 
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Onsager reciprocal relations now imply that the mag- 
netization dynamics induced by charge and spin cur- 
rents is given by Mj, = -[L (-m)^^]^fc(V^ - V^) + 
[L(— m)^|^^]cK/c/ifi^^/fi. Summarized, the generalized 
Landau-Lifshitz- Gilbert equation is 



m -7m X Htot + m X _ _ /_ m , (5) 

jMsV 

where 7 is (minus) the gyromagnetic ratio, Htot = Heff + 
+ Hf + Hc, and 

H, (m) ^ Tr[(Q^ - Q^) (-m)] (F^ - F^) /2, (6) 

Hf (m) = -Tt[& • Axf ^^^Q^^^) (-m)]/2e . (7) 

Heff = —dmF[^]/VMs, and we have introduced the 
Gilbert damping tensor in terms of the scattering matrix: 
[G(m)],, = 72?.J?e{Tr[(5^ jt)]}/4^ ^11^. Eqs. (§, 

(|6|, and ([7| are our main results and capture general 
charge-current induced magnetization-torque effects in- 
duced by any spin-orbit coupling, interface torques as 
well as bulk torque effects. The conventional spin- 
transfer torques can only be included with normal leads, 
where the spin is conserved. 

To check our formalism, let us now consider a system 
without spin-flip processes and derive the spin transfer 
torque on m from a spin accumulation fi^. For such a 
system, we expect our theory to be equivalent to mag- 
netoelectronic circuit theory [4]. In absence of spin-orbit 
interactions the scattering matrix is written as Smn,ii' = 

^InnM'^^^^tnM'^^ [S; , whcrc = (ao±(T-m)/2 is the 
matrix that projects out the spinor component for spin 
up (down) along the magnetization direction m. Substi- 
tuting this scattering matrix expression into Eq. ([5| gives 
the torque: 
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[gl^m X /L^f X m + gj^/j.^ x m ) . (8) 
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and assumed the ferromagnet is wider 
than the transverse spin coherence length, implying 
Emntlnn,LRittnn,LRy ^ ^' Eq. (§ exactly agrees with 
magnetoelectronic circuit theory [J, but is now derived 
from spin pumping via the Onsager relations. 

We will in the following apply our theory to a dis- 
ordered GaAs|(Ga,Mn)As|GaAs system, and investi- 
gate the magnetization dynamics in the ferromagnetic 
(Ga,Mn)As layer induced by an unpolarized charge cur- 
rent. We assume that the free-energy density depends on 
m as: 

T[m] — Kci {mlrriy + rnlml + rriyml) + Kufnl . (9) 

Here, Kd is the lowest-order cubic anisotropy con- 
stant [12]. We assume the system is grown on a 



GaAs substrate that induces compressive strain in the 
(Ga,Mn)As layer. Strain breaks the cubic symmetry. 
For growth direction [001] this results in the unaxial 
anisotropy energy Kuinl We use the anisotropy 

constants 2Kci/Ms 0.12 T and 2Ku/Ms 0.35 T 0. 

To model the band structure of (Ga,Mn)As we use the 
following Hamiltonian [ 12] : 



H = Hl^ h(r) • J + V{t) + H,, 



(10) 



Here, Hl is the 4x4 Luttinger Hamiltonian for 
zincblende semiconductors in the spherical approxima- 
tion. J is a vector of 4 x 4 spin matrices for J = 3/2 
spins. The h • J term is a mean field approximation 
of the exchange interaction between the itinerant holes 
and the local magnetic moment of the Mn dopants. The 
exchange field h is antiparallel to the magnetization di- 
rection m. V{t) = Vid{T — Hi) is the impurity po- 
tential, where is the position of impurity z, and Vi 
its scattering strength, which is randomly and uniformly 
distributed in the interval [— ^0/2,^/2]. Compressive 
strain in the (Ga,Mn)As layer induces the term i^strain- 
For growth direction [001], it is given by (i^strain)^^ = 
-6eax {J', - JV3)_ - ^C, {Jydy - J,5,)_ [l3], where 
b = -1.7 eV [13 , 6ax = (ezz - 6xx) = -0.005 El [l3], 
Cn^n are diagonal elements of the stress tensor, C4 = 
5 eV A [13], and 

Smn is the Kronecker delta. The GaAs 
leads are model by Hl^ only. 

We calculate the scattering matrix numerically us- 
ing a stabilized transfer matrix method. We consider 
a discrete (Ga,Mn)As layer with transverse dimensions 
Lx = 19 nm, Ly = 50 nm, and = 15 nm. The lat- 
tice constant is 1 nm, much less than the Fermi wave- 
length Ai? ~ 10 nm. The Fermi energy is 77 meV when 
measured from the lowest subband in (Ga,Mn)As. To 
estimate a typical saturation value of the magnetization 
we use Ms = lOl^lkx / a q^^j^^ p^, with x = 0.05 as the 
doping level, and where acaAs is the lattice constant 
for GaAs. The Luttinger parameters are 71 = 7.0 and 
72 = 2.5, |h| = 0.032 eV, and Vq = 0.4 eV for a dif- 
fusive system. Further details about the model and the 
numerical method can be found in Refs. [12 . 

We solve Eq. ([5| for the magnetization direction m 
in (Ga,Mn)As with Hc(m) and G(m) computed numer- 
ically from the scattering matrix. The scattering ma- 
trix is, for each magnetization direction, calculated for 
the (Ga,Mn)As region plus one lattice point into each of 
the GaAs layers. Heff is found from the free energy in 
Eq. ([9| . The magnetization starts initially along the easy 
X axis when an unpolarized current is applied to the sys- 
tem along the y axis. We find for a diffusive system with 
mean free path / ~ 25 nm that current densities as low 
as 5 X 10^ A/cm^ give magnetization switching between 
the X and y easy axis. 

cx m X 



Fig. |2|a) shows the vector field (m) 
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quantitative agreement with our results. In comparison, 
in metallic spin valve structures one typically observes 
critical spin currents of the order 10^ — 10^ A/cm^ [T]. 
Thus, charge current torque is not a negligible effect com- 
pared to the conventional spin transfer torque, and will 
play an important role in modelling spintronic devices 
made out of ferromagnetic semiconductors and other sys- 
tems with strong spin-orbit coupling. 

In conclusion, we have developed a general scattering 
matrix theory for magnetization dynamics that treats 
spin transfer torque and charge current torque induced by 
spin-orbit interaction. We apply our theory to a layered 
GaAs|(Ga,Mn)As|GaAs system and show that unpolar- 
ized charge currents switches the magnetization direction 
for current densities as low as 5 x 10^ A/cm^, m agree- 
ment with recent experiments. 

We thank Anh Kiet Nguyen for developing the numer- 
ical transfer matrix code. This work was supported in 
part by computing time through the Notur project and 
EC Contract IST-033749 "DynaMax". 
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FIG. 2: (Color online) (a) Vector field plot of (m 
m X He for a diffusive system with mean free path / ^ 
25 nm. Here, 6 and (f) are respectively the polar and az- 
imuth angles describing the local magnetization direction 
m = (sin^cos(/), sin^sin(/), cos 6). (b) Solution of Eq. ^ 
when applying a current density of 

y- 



-1 X 10 A/ CUT along 



azimuth angles and (j). The strong circular flow around 
the y axis, shows that He acts as an effective field along 
y. This torque-field arises solely from the strain term 
of the Hamiltonian in Eq. (10) [SlfMlIT^. In presence 



of an electric field along iYstrain oc Jyky {kyi Crys- 
tal momentum along y) tends to polarize the spin 3/2 
holes along the y axis. In combination with scattering at 
the GaAs|(Ga,Mn)As interfaces and impurities, this in- 
duces an out-of-equilibrium spin-density along y, giving 
the torque field shown in Fig. |2|a). 

Fig. |2jb) shows a solution of Eq. ([5| with the magne- 
tization starting initially along x when a current density 
of —1 X 10^ A/cm^ is applied at time t = along y. The 
evolution of m can be understood as follows: The charge 
field He rotates m a small angel around the y axis, giving 
a rriz component to the magnetization. This induces an 
effective field component along arising from the large 
unaxial anisotropy energy KuTnl in Eq. ([o]), that rotates 
m toward the y easy axis where the magnetization equi- 
librates. The inverse process is obtained by applying a 
current along x. 

Very recently, Ref. [3 reported an experimental obser- 
vation of charge current induced magnetization switching 
in (Ga,Mn)As for current densities just below 10^ A/cm^, 
while a theoretical study in Ref. calculated critical 
currents of the order 10^ — 10^ A/cm^. This is in semi- 
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